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ABSTRACT 


° 


The United States Coast Guard 1s organized by Atlantic and 
Pacific areas, which are further subdivided into districts. 
Each district assigns cutters (ships) of length 180 feet or 
less into weekly statuses. The resulting cutter schedules 
reflect the district’s level of readiness to respond to such 
emergencies as search and rescue, law enforcement, and 
pollution response. The First Coast Guard District has one of 
the largest scheduling problems, assigning each of 16 cutters 
to one of six weekly Statuses. The First District’s quarterly 
schedules must adhere to a number of guidelines which ensure 
patrol coverage, enforce equitable distribution of patrols, 
and restrict consecutive cutter statuses. This thesis 
formulates and solves the quarterly scheduling problem as an 
elastic mixed integer linear program. Face valid schedules, 
which are superior to actual schedules for all measures of 
effectiveness considered, are obtained within 15 minutes ona 
486/33 Mhz personal computer using a commercially available 
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I. INTRODUCTION 


A. COAST GUARD DISTRICT CUTTERS 

United States Coast Guard (USCG) District Commands 
schedule cutters of length 180 feet or less into weekly 
statuses. The First Coast Guard District has one of the 
largest scheduling problems, assigning each of 16 cutters to 
one of six statuses. The First District’s quarterly schedules 
must adhere to a number of guidelines which ensure patrol 
coverage, enforce equitable distribution of patrols, and 
restrict consecutive cutter statuses. 

Each cutter is scheduled by week to be either in an 
operational or maintenance status. The operational statuses 
are referred to as "bravo statuses". The most common bravo 
statuses used by the districts are: 


@® Bravo-24 (B-24) - The cutter must be underway within 24 
hours of notification to sail, 


@® Bravo-12 (B-12) - Within 12 hours, 
@® Bravo-6 (B-6) - Six hours, 
@® Bravo-2 (B-2) - Two hours. 


Cutters are also assigned various maintenance statuses, of 
which the most commonly used is "C" for "charlie periods". 

Bravo-2 and Bravo-6 statuses are used by the district for 
cutters which are designated as the first response vessels, or 


cutters which are required to patrol a certain area. Bravo-24 


and Bravo-12 statuses are assigned to allow a cutter to 
achieve routine maintenance or training. A charlie status is 
assigned to allow for more in-depth maintenance on the 
mechanical plant of the cutter and indicates the vessel will 
not be called upon to respond to anything but the gravest 
emergencies. 

In addition to the operational statuses, the cutter 
schedules also include yard and dockside availabilities. 
These availabilities are maintenance periods of longer 
duration than charlie periods, and are generally planned in 
advance through the USCG’s Maintenance and Logistics Command. 

Since Bravo-2 and Bravo-6 statuses fatigue both the crew 
and the cutter, the district scheduler attempts to rotate 
these statuses evenly among all available cutters. The 
scheduler also attempts to interrupt these statuses with less 
demanding statuses such as Bravo-12 and Bravo-24. The 
scheduler’s job is further complicated by the fact each class 
OE cutter SiS) required so haveu. certain number of weeks 
dedicated for maintenance, either charlie periods or dockside 
availabilities. There are also other demands on a cutter’s 
time, such as public appearances and patrolling special 


events. 


B. CUTTER SCHEDULES AND THE FIRST DISTRICT 
We define a "Schedule" as a specification of the weekly 


Status for all 16 cutters, which includes a specific patrol 


area with a patrol status. A quarterly (yearly) schedule is 
defined as a schedule for 13 (52) consecutive weeks. A 
"rough" schedule is defined as a schedule that does not 
indicate a specific patrol area with a patrol status. 

1. First District Cutters and Patrols 

The USCG’s First District stretches from Maine to 
halfway down New Jersey, and the district boundaries extend 
seaward to cover the ocean areas claimed by the United States. 
The district office is located in Boston, Massachusetts. 

The district scheduler is responsible for scheduling 
16 cutters from the following classes: 

@® 140 WIGB (icebreaking tug), 

@ 110 WPB (patrol boat), 

® 82 WPB (patrol boat). 
These cutters can be scheduled to four different patrol 
statuses. Each patrol status represents a different 
geographical area in the district which are designated: 
B-2NY, B-2G, B-2M, and B-2SAR. Not all of the cutter classes 
are eligible for every type of patrol year round. 

The 82 foot patrol boats cannot be deployed for search 
and rescue (SAR) during the winter months. This is because 
severe winter weather and icing conditions can cause stability 
problems for the cutters. The 140 foot WIGB cutters are also 
not available during the winter for patrols, because they are 


used to clear navigable harbors and rivers of ice. 


2. Current Manual Scheduling Practices 

Currently the district scheduler creates a rough 
annual schedule in September with paper and pencil. The 
scheduler then inputs the rough schedule to a spreadsheet. 
This process takes approximately two to three working days, 
and only indicates if a cutter will have a patrol status. 

Using this as a guide, the scheduler’s office 
subsequently assigns patrol areas (B-2NY, B-2G, B-2M, B-2SAR) 
to the designated patrol cutters, usually on a monthly basis. 
The job of creating an annual schedule has heretofore been 
considered too difficult for the scheduler to attempt. 

The rough annual schedule is constantly changing, 
because ships suffer casualties (break down) or requirements 
change. This makes the rough schedule developed in September 
less applicable the further into the future it is projected. 
Towards the end of each quarter, or more frequently should the 
need arise, the First District generates a new rough schedule 
for the remaining quarters. This is a very labor intensive 
practice with no clear measure of whether one rough schedule 
1s better than another. Essentially, if the scheduler and his 
assistant can fill the majority of the required patrols, the 


rough schedule is approved. 


C. OBJECTIVE OF CURRENT RESEARCH 
The objective of this thesis is to develop an 


optimization-based model to produce an annual schedule for the 


First Coast Guard District. Since the scheduler evaluates 
schedules quarterly, the problem is formulated on a quarterly 
basis. This allows the scheduler to generate an annual 
schedule, yet the individual quarterly schedules can easily be 
updated. 

Another goal of this thesis is to formulate and solve the 
model using GAMS (Brooke et al. (1988)) and XA (Byer (1987)). 
The use of commercially available software and a standard 
personal computer will allow the Coast Guard to quickly 
implement the model for a reasonable cost. Further, each 
Coast Guard District has its own constraints for scheduling 
cutters. Using commercially available software allows better 
support for the Coast Guard, and rapid adaptation of the basic 
model to the unique characteristics of the individual 


district’s own problems. 


D. THESIS OUTLINE 

Chapter II surveys related ship scheduling, models. In 
Chapter III, a mixed integer linear program is developed with 
detailed discussion of the costs, penalties and measures of 
effectiveness used. Using the model of Chapter III, 
computational performance for a year’s worth of actual data 
from the First District is presented in Chapter IV. 
Conclusions are provided in Chapter V. Appendix A contains a 
face valid annual schedule produced by the model, and Appendix 


B provides the GAMS program used to generate the model. 


Bae PREVIOUS SHIP SCHEDULING RESEARCH 


A. GENERAL SHIP SCHEDULING MODELS 

Most ship scheduling models reported in the operations 
research literature address problems faced by commercial 
shipping companies. Ronen (1983) provides an excellent review 
of the models which have been proposed. He discusses the 
variety and complexity of ship scheduling problems and 
proposes a model classification scheme. 

The models Ronen discusses are essentially concerned with 
a fleet of ships moving goods from one or more supply points 
to various demand points where the objective either minimizes 
the number of ships required in the fleet or minimizes 
transportation costs. It is clear the objectives and 
constraints used for commercial shipping companies are not 
directly applicable to the problem faced by the USCG 


districts. 


B. MILITARY SHIP SCHEDULING MODELS 
1. SeaLift Models 
Few military ship scheduling models have been 
developed. Of these, the most commonly seen formulations are 
Similar to the problems addressed for commercial purposes. 
Lally (1987) and Lima (1988) formulate problems for the 


Emergency Deployment Agency. The main objective of these 


models is to move as much material from several embarkation 
ports to disembarkation ports in as little time as possible. 
2. U.S. Navy Atlantic Fleet Scheduling Problem 

Ratliff (1981) first explored the possibilities of 
using an integer program for scheduling the Navy’s Atlantic 
Fleet. This work was later extended in Ratliff and Nulty 
(1986) where a network-based model is developed for scheduling 
the Navy’s Atlantic Fleet. In Ratliff and Nulty (1986) each 
individual ship’s schedule is viewed as a network and solved 
aS a longest path problem. These longest path problems are 
solved repeatedly, changing certain parameters, until all of 
the specified requirements are satisfied. Computational 
results are reported that indicate sample problems can be 
solved to near optimal with between five to twenty iterations 
of the algorithm; however, no indications of the computation 
requirements are provided. 

Goodman (1985), followed by Brown et al. (1990), 
develop an extremely efficient algorithm for scheduling 
surface combatants of the Atlantic Fleet, which they entitle 
CPSKED. The CPSKED model uses a column generation technique 
that includes all possible feasible schedules. It then uses 
an elastic set partitioning model to select the best set of 
candidate schedules. The authors report solutions to the full 
scale 111-ship annual scheduling problem in under 2 minutes on 


an IBM 3033 AP computer. While conceptually the problem of 


scheduling combatants for the Atlantic Fleet is similar to 
that faced by USCG cutters, CPSKED is concerned with matching 
ship capabilities (armament, communication) with missions of 
varied durations. These considerations are not applicable to 
this thesis. 
3. Quadratic Assignment/Linear Programming Model 

Sibre (1977) developed the only model found in the 
literature which specifically addresses a Coast Guard 
scheduling problem. Sibre’s Naval Postgraduate School thesis 
contains a mathematical model to schedule the Pacific Area’s 
Hamilton Class High Endurance Cutters. The model is solved 
using a variation of the Quadratic Assignment/Linear 
Programming model developed by Geoffrion and Graves (1976). 
The Quadratic Assignment model contains the guidelines for the 
number of cutters needed in a patrol area, the cutters 
available for a specific assignment, the transition cost of 
cutters shifting from one type of assignment to another, and 
the costs for violating starting and completion dates. The 
linear programming model is used primarily to determine the 
length of the patrols. The model develops a schedule and then 
uses a heuristic approach to alter the schedule to ascertain 
1f any improvement can be achieved. 

Sibre’s model cannot be easily adapted to the district 
cutter scheduling problem because of fundamental differences 


between the Area scheduling problem and the district 


scheduling problem. The Area scheduling problem has patrol 
lengths of variable duration. fiMeseCont- rast, sene Gistrict 
schedules patrols of exactly one week. A major concern of 
Sibre’s model is the extensive transit times faced by the High 
Endurance Cutters between patrol areas. Since the High 
Endurance Cutters transit to patrol areas far away from their 
homeport, days away from homeport became a good measure of 
effectiveness. In the problem under consideration in this 
thesis, neither of these factors are applicable. 
4. Area Scheduling Model In Development 

Currently the Coast Guard’s Research and Development 
Center in Groton, Connecticut, has a project to develop an 
optimization based model to create annual schedules for the 
Atlantic Area. The R&D Center contracted a group from Brunel 
University in London, England, led by Dr. G. Mitra to 
formulate and solve the problem. The problem is to schedule 
Coast Guard Medium and High Endurance Cutters to various 
assignments on the East Coast of the United States. The 
problem calls for daily time resolution in the schedule 
produced, accounting for transit time, and certain constraints 
imposed by the Atlantic Area Command. While this research has 
not been published (see Darby-Dowman et al. (1992)), the 
problem is similar to that addressed by Brown et al. (1990) 
and Sibre (1977) and therefore possessed many characteristics 


which make it not directly applicable to this thesis. 


BA a Be FIRST DISTRICT SCHEDULING MODEL 


A. PROBLEM DESCRIPTION AND FORMULATION 

The problem of creating a quarterly schedule for the First 
Coast Guard District is formulated as an elastic mixed integer 
linear program. The formulation ensures all patrol 
requirements are satisfied. The elastic portion of the model 
ensures that requirements are minimally violated when they 
cannot be feasibly satisfied. The model assigns 16 cutters to 
one of the following statuses every week of the quarter: 


@® B-2M - The cutter must be within the geographic boundaries 
of the northern patrol area designated by the letter M, 


@® B-2G - The cutter must be within the central patrol area 


G, 

@® B-2NY - The cutter must be within the southern patrol area 
NY, 

@® B-2SAR - A patrol status indicating the cutter is a 


designated search and rescue vessel, 

@® B-12 - The First District uses this status for any week a 
cutter is not scheduled for a specific patrol or 
maintenance period, 


ec - This status indicates a cutter is assigned a charlie 
period, 


@® UNAVAL - This status is assigned when a cutter is not 
available for patrol assignments. 


The model allows cutters to be manually specified as 
unavailable and makes sure the required number of "charlie 


periods" for each cutter are assigned. Additionally, the 
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model enforces all of the First District’s policies regarding 


cutter scheduling, which include: 


Charlie periods should be in two-week blocks, 


Cutter’s cannot have more than two consecutive patrol 
statuses, 


Cutter’s cannot be assigned consecutive SAR statuses, 


82 WPBs cannot be assigned SAR patrols during the Fall and 
Winter Quarters, 


WITGB class cannot be assigned SAR patrols, 


Patrols should be assigned to cutters equitably. 


CutS (CUTTER SCHEDULER) 


The formulation of CutS (Cutter Scheduler) is presented 


below after the introduction of appropriate notation. 





Indices: 
ies Cutter; 
k = statuses (B-2M, B-2G, B-2NY, B-2SAR, C); 
t = week the cutter assumes the patrol status. 
Data 
COST, = COsE OF scheduling CUmGer:1 FO patrol; 
aL eige eclobljomaat is available for patrol 
SHIPAVAL,, = during week t; 
QO otherwise; 
REQ, = required number of cutters for status k 


excluding <¢™- 


LL: 


CHARLIE, 


RPEN, 


CPEN 


FPEN 


FAIRLO, 
FAIRHI, 


VARIABLES: 


es 


Minimum number of weeks of charlie status 
required for cutter 1; 

penalty for not meeting the required 
number of cutters (REQ,) for patrol status 
k; 

penalty for violating charlie status 
constraints; 

penalty EOr violating fairness 
constraints; 

minimum number of patrols for cutter 1; 


maximum number of patrols for cutter i. 


if cutter iis assigned status k in week t; 


otherwise; 


elastic variable measuring deviation from 


required patrol statuses; 


elastic variable for violating the consecutive 


charlie period requirement; 


elastic variable for violating lLimitWonmne 


more than two consecutive charlie periods; 


elastic variable for violating fairness 


constraints. 


LZ 


EQUATIONS : 


Minimize }/ p> Yo (XineXCOST,;, + DREO,,xRPEN, + DCON,,xCPEN + 
oi c 


D2C,,.xCPEN + DF,xFPEN) 


Subject to: 


(1) > Xize = REO, - DREO,, Vi) (ie ac, 6) 


(2) ee ee (7) 


(3) \> Xice 2 CHARLIE, va 
C 
(5) Ke ee 2D eC eV (tec rg) 
(6) DO ae EATREO, DF, VY i 
kee Ct 
(7) Yd Xie, S FAIRHI, + DF; PS 
kec C 
(8) eee eee Mae + XG.) Ss 270 Va ee) 
kee 
(9) Kip-asare * “ip-2sare-1 S 1.0 vo (2, €> 1) 


CONSTRAINT EXPLANATION: 


aly, A minimum number of cutters must be assigned to 


patrol status k. 


gb 


(2) Each cutter cannot be assigned more than one status 
for each week. Any cutter not assigned a status is placed in 
B-12. 

(3) Each cutter must have at least a minimum number of 
required charlie periods. 

(4) Charlie periods must be consecutive. 


(5) Consecutive charlie periods should not exceed two 


(6) Each cutter must have a minimum number of patrols. 

(7) Each cutter cannot exceed a maximum number of 
patrols. 

(8) A cutter cannot be assigned more than two consecutive 
patrol statuses. 

(9) A cutter cannot be assigned consecutive B-2SAR 


statuses. 


C. CONSTRAINT DESCRIPTION 

Constraints (4) and (5) require charlie periods to be 
assigned in two-week blocks. The First District Scheduler has 
not effectively been able to manually assign charlie periods 
in this manner due to the complexity of the task. 

One of the goals of CutS is to create quarterly 
schedules which distribute the number of patrols between the 
cutters equitably; constraints (7) and (8) are designed to 
achieve this goal. The parameters FAIRHI, -and FAIRLO, are 


calculated based on a number, FAIR, which is the total number 
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pereatired patrols divided by the number of cutters. FAIRHI, 
1s established by simply adding two to the number FAIR. 
Allowing the number of patrols assigned to each cutter to 
deviate by two empirically gives the model flexibility, while 
Still adhering to the First District Scheduler’s requirement 
to equitably distribute the patrol statuses. FAIRLO, is the 
Minimum of FAIR and the maximum number of patrol assignments 
possible for each cutter. The maximum is easily determined by 
summing the number of weeks the cutter is available and 
subtracting the number of required charlie periods. 

The last two constraints, (8) and (9), express First 
District policy, requiring that no cutter have more than two 
consecutive patrol periods. The District’s policy also 


dictates cutters should not have consecutive SAR statuses. 


D. COSTS, PENALTIES, AND MEASURES OF EFFECTIVENESS 
The difficulty establishing meaningful measures of 

effectiveness for ship scheduling problems is well documented 
by Soland (1982). The measures of effectiveness used for CutS 
are essentially the same as those employed by the First 
District Scheduler: 

@® Minimize the transit time to cutter’s patrol areas, 

@® Minimize the number of required patrol statuses missed, 

® Equitably distribute the patrol statuses among cutters. 
In addition to the above measures of effectiveness, the amount 


of time necessary to create quarterly and annual schedules is 
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considered in the computational results reported in Chapter 
IVR 

The costs and penalties for CutS stem from the measures of 
effectiveness and from discussions with the First District 
Scheduler. The cost of assigning a cutter to a patrol is the 
transit time of the cutter to the patrol area. Accordingly, 
the transit time matrix shown in Table 1 was developed. The 
entries in Table 1 indicate the hours required to transit from 
a cutter’s homeport to the patrol area at the cutter’s normal 
cruising speed. For example, it takes the ADAK 17 hours to 
transit from its homeport to the B-2M patrol area. Even 
though cutters generally remain in their homeport for their 
charlie periods, a transit time of two is used for any cutter 
assigned to "C". The number two is used because it is less 
than the smallest transit time in Table 1, but it is greater 
than the zero used for the B-12 status. This relationship 
ensures that only the minimum number of charlie periods are 
assigned to each cutter. 

CutS uses the following penalties, which are derived from 
the transit time matrix: 

@® 40 for missing a B-2M, B-2NY, or B-2SAR status, 
@e 30 for missing a B-2G status, 


@® 25 for each unit of violation from the fairness 
constraints, 


@® 40 for violating charlie requirements. 


LG 


The scheduler desires all patrol requirements to be satisfied, 
regardless of the transit time for a cutter. This results in 
the penalty of 40 for missing a required status (a value 
slightly greater than any of the values in Table 1). Failing 
to fill a B-2G status is penalized less than the other patrol 
Statuses, because two cutters are normally assigned to this 
patrol area. If the scheduler is unable to meet all of the 
required patrol statuses, the first status to be unfilled will 
likely be one of the B-2G patrols. A penalty of 25 is used 
for violating the fairness constraints, since the scheduler 
indicates it is better to assign a cutter extra patrols than 
to leave an area uncovered. However, the scheduler would not 
assign an extra patrol to a cutter which is very far away from 
the cutter’s homeport. The penalty of 25 allows CutS numerous 
alternatives to reasonably violate the fairness constraint 


prior to exceeding the penalty for missing a B-2G status. 


ey 


TABLE 1 
TRANSIT TIME MATRIX 


The values in this table represent the transit time in hours 
for a cutter to reach a patrol area. These values are used in 
the objective function of sGiles: CutS’ objective function 
minimizes the total quarterly transit time used by cutters to 
reach patrol areas. 
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Iv. COMPUTATIONAL EXPERIENCE 


A. TEST PROBLEMS 

The First District Scheduler provided the District’s rough 
quarterly schedules for the last two quarters of fiscal year 
1991 and the first two quarters of fiscal year 1992. The 
rough quarterly schedules yielded the following information, 
which 1S summarized in Table 2: 


@® the weeks cutters were not available to be assigned any 
statuses, 


@® the number of charlie periods each cutter was assigned, 


@® the number of required patrol statuses which were missed 
during the quarter, _— 


@® the number of patrol statuses assigned to each cutter. 
The First District also provided a 12 week schedule for the 
second quarter of fiscal year 1992. 

The quarterly versions of CutS were run with the 
information shown in the first two columns of Table 3. While 
the total number of charlie periods per quarter were virtually 
the same between CutS and the provided rough quarterly 
schedules, some of the individual cutter’s charlie periods 
were modified by the addition or subtraction of one period. 
These modifications were conducted to establish an even 
charlie period requirement for each cutter and thereby more 


applicably model constraints (4) and (5). 
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TABLE 2 
SUMMARIZED INFORMATION FROM FIRST DISTRICT’S SCHEDULES 


"Weeks of unavailable" represent the total number of weeks 
cutters were not available to receive patrol assignments. 
"Weeks of charlie" is the number of charlie periods given to 
the cutters during the quarter. "Missed statuses" is the 
number of required patrols the manually created schedule was 
unable to fill. 


| Weeks of Weeks of Missed 
| Unavailable Charlie Statuses 
Sealy, Suastas fe} og 


Spring Quarter a Yd 
[Summer quarter | 53 | 63 | 







TABLE 3 
SUMMARIZED TEST PROBLEM DATA INPUT 


Using essentially the same initial conditions as Table 2, CutS 
had significantly fewer missed statuses. 


Weeks of Weeks of Missed 
———| a = —" 


Fall | Fall Quarter _| 








|/Winter Quarter 


spring quarcest [5 _46. gualeemeaeeen 
[Summer Quarter | 53 | 62 | 
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All computational results reported in this thesis were 
obtained using a 486/33 Mhz personal computer and the 


commercial linear/integer programming solver XA. 


B. SOLUTION TIMES 
A basic measure of effectiveness is the time required to 
obtain feasible quarterly schedules. The schedules generated 
each quarter using the information of Table 3 are included in 
Appendix A. These quarterly schedules were shown to the First 
District Scheduler for his critique. The scheduler and his 
assistant reported the schedules are feasible and of good 
quality. 
1. Time Necessary to Develop Schedules 
The scheduler and his assistant take approximately two 
to three working days to develop a rough annual schedule. The 
amount of time CutS requires to develop each quarterly 
schedule within specified tolerances of optimality are shown 
in Table 4. Generation of an annual schedule within 5% of 
optimal by running the quarterly versions of CutS 
consecutively, including the required inputs for each model, 
can easily be accomplished in under two hours. This is a vast 
improvement over the two to three days required by the 
district scheduler. 
The Fall and Winter versions of CutS consistently take 
longer to solve than the Spring and Summer versions of the 


model. This is partly because of the restrictions imposed on 
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using the WTGB and 82 WPB class of cutters. It is also the 
result of the information provided by the First District, 
which was used as the initial inputs for CutS. The number of 
free weeks cutters had available for patrols during the Fall 
and Winter are Significantly less than the weeks available for 
the Spring and Summer, as can be seen in Table 2. These two 
factors make the Fall and Winter models have comparatively 
less scheduling flexibility than the Spring and Summer. 

The solution times required to guarantee an optimal 
solution, or a solution within 1% of optimal, dramatically 
increase for more restrictive quarters as compared to solving 
within 5% of optimal. However, the quality of the quarterly 
schedules, as indicated by the objective function values shown 
in Table 5, do not display dramatic improvement for the test 
problem considered. For example, the time necessary for the 
Fall version of CutS to create a quarterly schedule guaranteed 
within 5% and 1% of optimal was 13.4 and 157.0 minutes 
respectively. The objective function value for the model, 
however, actually increased from 386 to 387. 

The model produces face valid quarterly schedules for 
anything within 10% of optimal. For solutions within 5% of 
optimal, CutS usually creates schedules which have not 
violated any of the constraints. The changes in the quarterly 
schedules above the 5% level are the result of CutS switching 
patrol senceaese which reduce the integer objective function 


value only slightly. 


22 


TABLE 4 
SOLUTION TIMES OF THE QUARTERLY MODELS 


The solution times are obtained using a 486/33 Mhz P.C. and 
show at any setting the ability of CutS to quickly obtain 
schedules which took more than two days to manually develop. 


Time in minutes to obtain solution within 
percentage from optimal 

ee ee ee ee 

ee ee eee se 

jowrse [5.2 [0.8 16.4 05.7 | 


SPRING 


la eee er 












TABLE 5 


OBJECTIVE FUNCTION VALUES 


This table demonstrates that objective function values 
Guaranteed to be between 5% and 1% of the optimal do not 
change appreciably, whereas Table 4 shows a significant 
increase in time. The 5% setting is therefore recommended. 
The linear program (LP) objective function value is an easily 
obtained bound of the best possible solution. 
















Objective function value for 
guaranteed within percentage from 


of LP 
optimal — Solution 
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2. Computational Difficulties Imposed by Constraints 

The requirement for charlie periods to be in two-week 
blocks (CutS constraints (4) and (5)) greatly increases the 
computational difficulty, and hence the amount of time 
necessary to obtain a solution. Because these constraints 
force CutS to assign charlie periods in two-week blocks, 
entering odd numbers for the minimum required charlie periods 
empirically makes it more difficult, if not impossible, for 
the model to solve within 10% of optimal. 

Test runs of CutS to investigate how much 
computational difficulty constraints (4) and (5) add are shown 
in Table 6, where a 14-hour time limit was imposed. Table 6 
compares computational performance with all cutters having 
requirements for an even number of charlie periods, 25% of the 
charlie requirements being odd numbers, and with constraints 
(4) and (5) removed. Table 6 dramatically portrays the 
computational difficulty imposed by these constraints and by 
entering an odd number of charlie periods. 

Note that requiring the total number of charlie 
periods to be an even number is not a eine assumption. 
The scheduler is given some flexibility when setting these 
requirements, and the total number is large enough that the 
addition or subtraction of one period is only a slight 


alteration. 
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C. UNFILLED PATROL STATUSES 
The most important measure of effectiveness is the number 
of required patrols missed. This criteria is essentially how 
the district scheduler compares the various rough quarterly 
schedules developed by hand; consequently, it is a logical 
TABLE 6 


CONSECUTIVE CHARLIE CONSTRAINTS AND COMPUTATIONAL 
REQUIREMENTS 


The computational difficulty imposed by the constraints on 
consecutive charlie periods is evident from the large increase 
in the solution times with the constraints added. The effect 
of using an odd number of required charlie periods is also 
demonstrated. 


Quarter Time in minutes to obtain golution within 10% 
ae | sc rom optimal or exceed time limit. 


All Even Number | 25% of Charlie | Constraints 
Charlie Periods | Periods Odd Removed 


Ce a | 840.0 


WINTER 840.0 


SPRING er an eT 
aa ee 


measure of how well CutS performs. The rough annual schedule 














developed by the First District leaves 19 statuses unfilled, 
Shown in Table 2. CutS annual schedule leaves only three 


Statuses unfilled, shown in Table 3, a marked improvement. 


D. EQUITABLE DISTRIBUTION OF PATROL STATUSES 
Another measure of effectiveness is how equitably the 
patrol statuses are distributed among the cutters. CutS is 


required to develop quarterly schedules which. fairly 
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distribute the number of patrol statuses assigned to each 
CUGEer (This is also a goal which the district scheduler 
uses.) The number of patrols assigned to each cutter per 
quarter, along with the yearly totals, are shown in Table 7 
for the rough quarterly schedules of the First Distrigel 
Table 8 displays the same information for the quarterly 
schedules generated by CutS. It is difficult to determine any 
Significant differences between the schedules by looking at 
the respective tables; consequently, the standard deviations 
are calculated. The standard deviation for yearly number of 
patrols assigned by the First District (the WTGB class is not 
included in the calculations) is 2.55 patrols, while the 
standard deviation for the CutS’ schedule is 2.40 patrols. A 
more impressive result is noted in the third and fourth 
quarters. These two quarters have the greatest number of 
cutters available for patrols, hence there is more flexibility 
in scheduling cutters to meet the requirements. The standard 
deviation for the First District’s third quarter schedule 
(Including the WIGB class) is 1.41 patrols, while that of CutS 
is .77 patrols. The fourth quarter’s schedules display 
Similar standard deviations with the district’s schedule being 
1.83 patrols, while the schedule proposed by CutS has a 
standard deviation of .57 patrols. CutS is able to generate 
schedules which more evenly distribute patrol statuses among 


the cutters. 
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TABLE 7 
NUMBER OF PATROLS MANUALLY ASSIGNED PER CUTTER, 
PERS oto eR LCr 


This table displays the number of patrols per quarter each 
cutter was assigned by the First District Scheduler. The 
table shows how equitably the manually developed schedule 
distributed the patrols by comparing entries within the same 
column. 


Po tse | ana | ara | ath | total 
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TABLE 8 
NUMBER OF PATROLS ASSIGNED PER CUTTER, OPTIMIZATION MODEL 


This table shows the number of patrols assigned to each cutter 
per quarter by CutS. This table indicates how equitably CutS 
was able to distribute the patrols by comparing entries within 
the same column. Any value below 4 was a result of limited 
cutter availability and should not be considered as 
inequitable. Contrasting these results to Table 7 shows CutS’ 
superiority to equitably distribute patrol assignments. 


a er ee 


mors. |e | = 1 7 nn 
pr-pontta | 4 | 4 | 4 | 4 | ae 


pr-rrancis | 6 | 6 | 4 | 4 | 20 
pr-gaccson | 4 | 5 | 4 | 3 | a6 
- HANN Js) 
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E. ASSIGNMENT OF PATROL STATUSES 

The final measure of effectiveness is the ability to 
assign patrols to cutters within the geographic area of,the 
cutter’s homeport. While the district scheduler makes patrol 
assignments with this in mind, the scheduler does not actually 
Calculate the transit times of the cutters. CutS, on the 
other hand, minimizes the objective function value, which is 
based on the total transit time of the cutters. The objective 
function value of the 12-week schedule produced by the First 
District for the Winter quarter is 458, while the objective 
function value of the schedule produced by CutS is 310. These 
objective function values are based on the cost of specific 
patrol assignments and on penalties ‘incurred for missing— 
required patrols. In order to get a better comparison between 
the First District and CutS, the four weeks for which no 
penalties are incurred for missing required patrols are 
compared. The objective function value associated with the 
First District’s schedule is 105, while that of the schedule 


produced by CutS is only 83. 


Zag 


Vi CONCLUSIONS 


CutS produces quarterly schedules in a reasonable amount 
of time on a personal computer which are superior to manually 
developed schedules in all areas of concern: schedules 
developed by CutS miss significantly fewer required patrols, 
the schedules assign patrols to cutters which are closer to 
the cutter’s homeport, and CutS’ schedules more equitably 
distribute the number of patrols each cutter receives. Cuts 
develops face-valid schedules which can be implemented without 
any changes; however, the real benefit of the model is that it 
will assist the district scheduler to quickly develop feasible 
schedules. This will give the scheduler an opportunity to 
produce quarterly schedules of a higher quality than has been 
previously possible. 

CutS was demonstrated for the First District Scheduler and 
members of his staff at the Coast Guard’s Research and 
Development Center. The scheduler brought a quarterly 
schedule he had spent many hours developing. He was not 
pleased with his results because of a high number of required 
patrol statuses missed. He was able to input ee data 
required for CutS and obtain a quarterly schedule, which did 
not miss any required patrols, in under 30 minutes. This 


demonstration resulted in a push for immediate implementation. 
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CutS can be improved by adding the capabilities of a 
spreadsheet to the model. The Coast Guard’s Research and 
Development Center has developed a spreadsheet and database to 
be used with the Area scheduling model under development. The 
Research and Development Center intends to adapt this 


interface for use with Cuts. 
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APPENDIX A 


QUARTERLY SCHEDULES DEVELOPED BY Cuts 


The following are quarterly schedules created by CutS, solved within 5% of 


optimal using the initial data listed in Table 3. 


The left hand margin 


has the cutter’s name, and the top margin lists the week the cutter is to 


have the assigned status. 


Week 1 commences on the Monday of the first 


week of the first quarter of a fiscal year and continues until the 
following Monday. Each quarterly schedule includes the last two weeks of 


the previous quarter. 


1ST QUARTER 


51 
ADAK B-2G 
WRANGEL UNAVAL 
SANIBEL UNAVAL 
MONOMOY Cc 
JBEFF-ISL C 
GRAND -ISL c 
BAIN-ISL B-2G 
PT-BONITA B-2NY 
PT-FRANCIS B-12 
PT - JACKSON S 
PT - HANNON B-2M 
PT-TURNER B-2SAR 
PT-WELLS UNAVAL 
PENOBSCOT UNAVAL 
STURGEON B-12 
THUNDER B-12 

4 
ADAK Cc 
WRANGEL B-2M 
SANIBEL B-2SAR 
MONOMOY UNAVAL 
JEFF-ISL c 
GRAND - ISL B-2G 
BAIN-ISL B-2NY 
PT-BONITA Cc 
PT-FRANCIS c 
PT -JACKSON B-2G 
PT -HANNON B-12 
PT -TURNER C 
PT -WELLS UNAVAL 
PENOBSCOT UNAVAL 
STURGEON UNAVAL 
THUNDER UNAVAL 
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i] 
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NO 
Q) 


g8 
BEB 
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ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JBEFF-ISL 
GRAND - ISL 
BAIN-ISL 
PT-BONITA 
PT-FRANCIS 
PT - JACKSON 


PENOBSCOT 
STURGEON 
THUNDER 


2ND QUARTER 


ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JEFF-ISL 
GRAND -ISL 
BAIN-ISL 
PT-BONITA 
PT-FRANCIS 
PT - JACKSON 
PT - HANNON 
PT - TURNER 
PT'-WELLS 
PENOBSCOT 
STURGEON 
THUNDER 


ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JEFF-ISL 
GRAND -ISL 
BAIN-ISL 
PT-BONITA 
PT-FRANCIS 
PT - JACKSON 
PT - HANNON 


PENOBSCOT 
STURGEON 
THUNDER 


B-2G 
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ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JBEFF-ISL 
GRAND-ISL 
BAIN-ISL 
PT-BONITA 
PT-FRANCIS 
PT - JACKSON 


PT-WELLS 
PENOBSCOT 
STURGEON 
THUNDER 


3RD QUARTER 


ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JEFF-ISL 
GRAND-ISL 
BAIN-ISL 
PT-BONITA 
PT- FRANCIS 
PT - JACKSON 
PT - HANNON 
PT -TURNER 
PT-WELLS 
PENOBSCOT 
STURGEON 
THUNDER 


ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JEFF-ISL 
GRAND-ISL 
BAIN-ISL 
PT-BONITA 
PT-FRANCIS 
PT - JACKSON 
PT - HANNON 
PT -TURNER 
PT-WELLS 
PENOBSCOT 
STURGEON 
THUNDER 


DANWWAAAAWWWAAN 


en 
pepe 


ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JEFF-ISL 
GRAND -ISL 
BAIN-ISL 
PT-BONITA 
PT-FRANCIS 
PT -JACKSON 
PT -HANNON 


PENOBSCOT 
STURGEON 
THUNDER 


4TH QUARTER 


ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JEFF-ISL 
GRAND -ISL 
BAIN-ISL 
PT-BONITA 
PT-FRANCIS 
PT - JACKSON 
PT - HANNON 
PT- TURNER 
PT-WELLS 
PENOBSCOT 
STURGEON 
THUNDER 


ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JEFF-ISL 
GRAND-ISL 
BAIN-ISL 
PT -BONITA 
PT-FRANCIS 
PT-JACKSON 
PT - HANNON 
PT-TURNER 
PT-WELLS 
PENOBSCOT 
STURGEON 
THUNDER 
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ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JBEFF-ISL 
GRAND -ISL 
BAIN-ISL 
PT-BONITA 
PT-FRANCIS 


PENOBSCOT 
STURGEON 
THUNDER 
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APPENDIX B 


FALL VERSION OF CutS 


STITLE FIRST DISTRICT SCHEDULING MODEL, 1st Quarter 
SOFFUPPER OFFSYMLIST OFFSYMXREF 


OPTIONS 
LIMCOL = 0, LIMROW = 0, SOLPRINT = OFF, RESLIM = 5000, 
OPTCR = 0.05, INTEGER1I = 1, ITERLIM = 100000 ; 


SONTEXT 
Original by: LT. Robert A. Farmer, USCG date 05/10/92 
Dr. Robert F. Dell, Naval Postgraduate 
School, Monterey, California 


Description: 

This is an optimization based model to solve the First District’s 
cutter scheduling problem. The program needs to be run with GAMS and an 
integer solver. All of the trial runs of this model were run with the 
integer solver XA. 

The indices greatly affect the number of variables created, which in 
turn affects the ability of the solver to find a reasonable solution. 
Care should be exercised when adding new elements to the indices. | 

This model will develop the first quarter schedule for the First Coast 
Guard District. The status of the cutters represents the following: 


B-2M - vessel assigned to the Gulf of Maine OPAREA 
B-2G - vessel assigned to the Georges Bank OPAREA 
B-2NY - vessel assigned to the New York Bight OPAREA 
B-2SAR - vessel assigned as the SAR standby cutter 

C ' + vessel in maintenance status 


**NOTE** B-12 is not included in the index for vessel statuses, but it 
will be included in the final output by the print statement at the end of 
this file. . 

Cutters are assigned by weeks, with week 1 being the first Monday in 
the first quarter of a fiscal year. No days-were included in the model 
for the transit time of the cutter to the OPAREA. 

If any changes are made to the sets listed below, it is imperative the 
same changes be made everywhere the set is used. For example if a new 
cutter the NEVERSAIL is added, NEVERSAIL must be added to SHIPAVAL, 
CHARLIE, etc. 

SOFFTEXT 


SETS I name of the cutter 
/ ADAK, WRANGEL, SANIBEL, MONOMOY, JEFF-ISL, GRAND-ISL, BAIN-ISL, 
PT-BONITA, PT-FRANCIS, PT-JACKSON, PT-HANNON, PT-TURNER, 
PT-WELLS, PENOBSCOT, STURGEON, THUNDER/ 


K status of cutter 
/B-2M, B-2G, B-2NY, B-2SAR, C/ 


ay 


T week the cutter assumes the assigned status 
fein Se) ake} 


SONTEXT 

The numbers in the cost table represent the transit time from 
the cutter’s homeport to the OPAREA of the patrol statuses. 
SOFFTEXT 


TABLE COST(I,k) cost of scheduling cutter I for patrol k 


B-2M B-2G B-2NY B-2SAR C 
ADAK 7 9 3 9 2 
WRANGEL 3 5 0 5 2 
SANIBEL 6 3 6 3 2 
MONOMOY 6 3 E 3 2 
JEFF -ISL 3 5 17 5 2 
GRAND -ISL 3 3 11 3 2 
BAIN-ISL 17 9 3 9 2 
PT-BONITA 18 7 4 7 2 2 
PT- FRANCIS 14 4 4 4 2 
PT - JACKSON 18 5 5 5 2 
PT - HANNON 4 ll 21 11 2 
PT - TURNER 1l 04 7 4 2 
PT-WELLS 14 5 4 5 2 
PENOBSCOT 3618 5 18 2 
STURGEON 36 «18 5 18 2 
THUNDER 6 18 36 18 2% 
SONTEXT 


The matrix SHIPAVAL indicates which weeks a cutter is 
available to be assigned patrols. This matrix needs to 
be updated to reflect the actual data prior to running the 
model. A 1 indicates the cutter is available for patrol 
assignment and a 0 indicates it is not. 

SOFFTEXT 


TABLE SHIPAVAL(I,T) weeks cutter 118 available for patrols 


S51 52 10. 11. Sl20 
ADAK 
WRANGEL 
SANIBEL 
MONOMOY 
JEFF-ISL 
GRAND - ISL 
BAIN-ISL 
PT-BONITA 
PT-FRANCIS 
PT - JACKSON 
PT - HANNON 
PT-TURNER 
PT-WELLS 
PENOBSCOT 
STURGEON 


THUNDER 


PHOOPRPHPHPRPRPRPRPRROORRK 
PROORPRPRPRPRPRFRRPRPROOR 

DOOORPKPRPRPRFPBRPRPOOOFRR 
DODOORPKRPRPRPRRPHPROOORN 
DOOORPKRPRPHRPRPRRRFPOOOOW 
DOOORPRPHBBPRPRRPRORREFR b& 
DOORPRPRPRPRPRPRPRPRRORREHUW 
DOOR PRP PRP HRPRBPORFRFRDA 
DOORPRPRPRPRPRPHPRPRHPORRRA 
DOORPHPRPPHPRPBPRPRPORRRO® 
DOOPKPRPPBPRPRPRPRPRPRRR Ww 
DOORPRPRPRPPRPRPRPRPRRRHO 
DOOPRPRPRPRPRPRPRPHPRPRPRPRR 
DOORPRPRPRPRPRPHPRPPRPHRRPRH 

DOORPRPRPRPRPRPRPRRPRPRRRW 


PARAMETER 


aid 


REQ(K) required number of patrol boats for each OPAREA 
/ B-2M 
B-2G 
B-2NY 
B-2SAR 
Cc 


OrRPPNEH 


/ 


STAT(K) patrol status for variables to clean up program 
/ B-2M 
B-2G 
B-2NY 
B-2SAR 
Cc 


OrFPRPH 


SONTEXT 

Penalties greatly affect the model. These penalties were 
developed from the specific cost matrix above. They were 
established to be reasonably larger than any of the costs 
of the patrols. The penalty for missing a B-2G status was 
slightly less than the other penalties because two cutters 
are assigned to this OPAREA. 


SOFFTEXT 
PENALTY (K) penalty for not meeting minimum requirements 
/ B-2M 40 
B-2G 30 


B-2NY 40 
B-2SAR 40 
cS 40 / 


SONTEXT 

The array below represents the number of weeks of charlie 
status a cutter is to be assigned during the quarter. Do 
not include any charlie statuses from previous quarters, nor 
any charlie statuses accounted for in the SHIPAVAL matrix. 
**NOTE** ALL NUMBERS IN THE ARRAY BELOW SHOULD BE EVEN. 

If more charlie periods are entered for the quarter than there 
are weeks available for patrol assignments as listed in 
SHIPAVAL, the model will not solve and report the problem 

is infeasible. 

SOFFTEXT 


CHARLIE (I) weeks of Charlie for the quarter for each cutter 
/ ADAK 


PT-BONITA 
PT-FRANCIS 
PT - JACKSON 
PT - HANNON 
PT - TURNER 
PT -WELLS 
PENOBSCOT 
STURGEON 
THUNDER . 


OOOKPNAVHNAAWAKAANAHN A 


3g 


WTGB (I) if a cutter is a WTGB 


PT-BONITA 
PT-FRANCIS 
PT - JACKSON 
PT - HANNON 
PT-TURNER 
PT -WELLS 
PENOBSCOT 
STURGEON 
THUNDER 


RPrRrRoOqodoOOooOoOoOoOOoOOoOOoOoOo Oo 


/ 


WPB(I) if a cutter is a 82 foot patrol boat 
/ 


PT - JACKSON 
PT - HANNON 
PT -TURNER 
= PT-WELLS 
PENOBSCOT 
STURGEON 
THUNDER 


6 
w 
O 
x 
4 
» 
COOFPFHFHFPHHOOOOOOCO 


SONTEXT 

These penalties are for the elastic variables used in 

the model. Once again penalties greatly affect the 

performance of the model. PEN2 was chosen to be comparable 

with the penalties for not filling all patrol statuses. 

PEN3 is the penalty for the elastic variable of the fairness constraints. 
It was chosen to allow a cutter to receive an extra patrol only if the 
OPAREA is reasonably close to the cutter’s homeport. 

SOFFTEXT : 


SCALARS PEN2 penalty for temporary charlie variables /40/ 
PEN3 penalty for fairness constraints /25S/ ; 


SONTEXT 
The following four parameters calculate the number of 
patrols each cutter should receive in order for there 
to be an equal workload. If for some reason a cutter 
should receive a different amount it can be easily done. 
For example, say the cutter WRANGEL was assigned more 
patrols last quarter than any other cutter, and it should be given rest 
this quarter. Right under 
FAIRHI(I) = FAIR + 2.0 ; 
Buta 


40 


FAIRHI ("WRANGEL") = maximum number of patrols desired ; 
**NOTE** If you assign a FAIRHI value which is less 
then the FAIRLO value, the model will not solve and 
report the problem is infeasible. To insure this 
does not happen, first run the model and find the value 
for FAIRLO from the LST file. Make the desired changes 
and run the model again. 
SOFFTEXT 


PARAMETER REQTOT total of the requirements ; 
REQTOT = SUM(K, REQ(K)) ; 


PARAMETER FAIR equal number of patrols for each cutter; 
FAIR = ROUND (REQTOT* (CARD(T)-2.0)/CARD(I)) ; 
PARAMETER FAIRLO(I) lower limit on number of patrols; 
FAIRLO(I) = SUM(T$(ORD(T) GT 2), SHIPAVAL(I,T)) - 
CHARLIE(I) ; 
PAIRLO(I)S(FAIRLO(I) LT 0) = 0 ; 
FPAIRLO(I)$(FAIRLO(I) GT FAIR) = FAIR ; 


PARAMETER FAIRHI(I) upper limit on number of patrols cutters can 
assigned ; 


FAIRHI(I) = FAIR + 2.0 ; 
VARIABLES 


X(I,K,T) 1 if cutter i assigned status k for week l 
otherwise 


me ir, T) elastic variable to allow for no more than 3 "C" in 
a row 

E(K,T) elastic variables for unfilled status 

TC2(I,T) elastic variable for consecutive charlie periods 
LIM (I) elastic variable for fairness constraints 


TOTCOST objective variable (total cost) ; 


BINARY VARIABLE X ; 
POSITIVE VARIABLE TC ; 
POSITIVE VARIABLE E ; 
POSITIVE VARIABLE TC2 ; 
POSITIVE VARIABLE LIM ; 


SONTEXT 

The .UP variables place an upper limit on the number of times 
the constraints in which the variable is used may be violated. 
SOFFTEXT 


me-UP (1, T) =. 226 , 
mer. UP(I) = 1.0 ; 


feet, K,"51") 
mae, K, "52") 


tou 


SONTEXT 

Below is the mechanism used to inform the model of the patrol 
assignments for the last two weeks of the previous quarter. 
Be sure the matrix SHIPAVAL also accurately reflects the last 
two weeks of the previous quarter. Only input the statuses 
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be 


listed under index k above. In particular "UNAVAL" and 

"B-12" should not be entered. A patrol can be fixed for any 

week of the quarter simply by using the same format as below 

with the cutter’s name, the patrol, and the week the cutter is 

to have the patrol. It is important to remember fixing variables can 
adversely affect the solution time of the model. 

**NOTE** The patrol assignments for the last two weeks of the previous 
quarter must be properly filled in along with SHIPAVAL in order for the 
model to return a reasonable schedule. 

SOFFTEXT 


X.FX ("ADAK", "B-2G", "51") 
X. FX ("MONOMOY "7 "Gi", "52") ; 
X,FA(“UEPE-25eee eC" "51" )) 8 oor; 
X.FX("“GRAND=TSES "CC", "S1") 2 220m 

. FX ("BAIN - ete ee = I-0 
("PT -BONITA", "B-2NY", "S1") 2 2 
("PT - JACKSON", tC See ae 
("PT-HANNON", "B-2M", "51") 
("PT-TURNER", "B-2SAR","S51" 
("ADAK", "B- 2G", "52") =~ 0; 
("MONOMOY", UBe2SAR 7752.) = 
(NUBE Batol" C" 1. a2) aaa On; 
("GRAND - ISL", CMs 152 ees 1 Ole 
("BAIN-ISL", ncn 152") sel.O ; 
("PpT- BONITA", TC", "S2") = ‘ 
("PT- FRANCIS", OB ZN 2 
(SPT -JACKSONI 4 C:, 92.) = 
("PT-HANNON", "B- 2M" "5 20%) : 
AS fotos "B-2SAR" /T) SWTGB (I) = 0.0 oo ———— 
X.UP(I,"B-2SAR", 7) $ (WEB (TZ) $ (ORD (T) GT .2))) = 020; 


uel 
“ie eH 


ehelehelelelelelelolelelelea 
SEE EEEEEEEEEEE 





EQUATIONS 


COVREQ (K, T) constraint to meet minimum coverage for status k 
MINCHARL (TI) constraint to meet minimum amount of charlie time 
WEEKLY (I, T) cannot assign a vessel more than 1 status 

CONCHARI1 (I,T) charlie periods need to be in consecutive periods 
CHARZ (1, 1} consecutive charlie periods should not exceed 2 
LOWLIM (T) sum of patrol status cannot be less than lower limit 
UPPERLIM (TI) sum of patrol status cannot be more than upper limit 
NOREPEAT(I,T) cannot have more than two consecutive patrols 

NOSAR (I, T) cannot have consecutive SAR patrols 

TOTCOSTE ; 


TOTCOSTS -.. 


TOTCOST =E= SUM((I,K,T)$((ORD(T) GT 2)$SHIPAVAL(I,T)), 
X(I,K,T) *COST(I,K)) + SUM((K,T)$(ORD(T) GT 2), 
PENALTY (K) *E(K,T)) + SUM(I, LIM(I) *PEN3) + 
SUM((I,T)$((ORD(T) GT 2)SSHIPAVAL(I,T)), PEN2*TC(I,T)) + 
SUM((I,T)$((ORD(T) GT 2)$SHIPAVAL(I,T)), TC2(I,T)*PEN2) ; 
COVREQ (K,T)$((ORD(T) GT 2) SSTAT(K)).. 
SUM(ISSHIPAVAL(I,T), X(I,K,T)) =E= REQ(K) - E(K,T) ; 


WEEKLY (I,T)$((ORD(T) GT 2)$SHIPAVAL(I,T)).. 
SUM(K, X(I,K,T)) =L= 1.0 ; 


MINCHARL (TI) . 


SUM (TS ( (ORD (T) GT 2) $SHIPAVAL(1I,T)), xX(1, "Gi 
=G= CHARLIE (i, 
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CONCHAR1 (I,T)$((ORD(T) GT 2)$(CHARLIE(I) GT 1) SSHIPAVAL(I,T)).. 


X(I,"C",T)SSHIPAVAL(I,T) - X(I,"C",T-1) SSHIPAVAL(I,T-1) 
X(I,"C",T+1) SSHIPAVAL(I,T+1) - TC2(1I,T) $SSHIPAVAL(I,T) 
=G= © ; 


CHAR2 (I,T)$((ORD(T) GT 2)$(CHARLIE(I) GT 1)$SHIPAVAL(I,T)).. 


X(I,"C",T)SSHIPAVAL(I,T) + X(I,"C",T-1)$SHIPAVAL(I,T-1) + 


X(I,"C",T-2) SSHIPAVAL(I,T-2) =L= 2 + TC(I,T) $SHIPAVAL(I,T) ; 


LOWLIM(TI) .. 


SUM((K,T)$((ORD(T) GT 2) $STAT(K) $SHIPAVAL(I,T)), X(I,K,T)) =G= 


FAIRLO(I) - LIM(I) ; 


UPPERLIM(TI).. 


SUM((K,T)$((ORD(T) GT 2)$STAT(K)$SHIPAVAL(I,T)), X(I,K,T)) =L= 


FAIRHI(I) + LIM(I) ; 


NOREPEAT (I,T)$((ORD(T) GT 2)$SHIPAVAL(I,T)).. 
SUM(KSSTAT(K), X(1I,K,T) $SHIPAVAL(I,T) + 
X(I,K,T-1) $SHIPAVAL(I,T-1) + X(1I,K,T-2) SSHIPAVAL(I,T-2) ) 


NOSAR(I,T)$((ORD(T) GT 2)$SHIPAVAL(I,T)$(WTGB(I) LT 1)).. 
X(I,"B-2SAR",T) S$SSHIPAVAL(I,T) + 
X(I,"B-2SAR",T-1)SSHIPAVAL(I,T-1) =L= 1 ; 


MODEL ONE /ALL/ ; 
SOLVE ONE USING MIP MINIMIZING TOTCOST ; 
DISPLAY X.L, E.L, TC.L, TC2.L, LIM.L, REQTOT, FAIRLO, FAIRHI, 


SONTEXT 

The parameters break down the total costs associated with 
the model into the coverage costs and the costs associated 
with penalties. 

SOFFTEXT 


PARAMETER COVCOST coverage cost of the objective function ; 
COVCOST = SUM((I,K,T), X.L(I,K,T) *COST(I,K)) ; 


PARAMETER PENCOST penalty cost of the objective function ; 
PENCOST = SUM((K,T), PENALTY(K) *E.L(K,T)) + SUM(I, 
LIM.L(I) *PEN3) + SUM((I,T), TC.L(I,T) *PEN2) + 
SUM((I,T), TC2.L(I,T) *PEN2) ; 


DISPLAY COVCOST, PENCOST ; 


SONTEXT 

The following code generates the quarter employment schedule 
report into the file SKED1D.put 

SOFFTEXT 


FILE SKEDID; 
PUT SKEDID; 
SKEDID.NW = 12; 
SKED1D.TW = 12; 
PUT ’1ST QUARTER’ // 

@15, LOOP(T, 

IF (ORD(T) LT 6, 
PUT TTL) : 
ie 
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azL= 2 


FAIR; 


LOOP(I, PUT / I.TL; 
LOOP (T, 
IF (ORD(T) LT 6, 
IF (SHIPAVAL(I,T) EQ 0, 
PUT ‘UNAVAL’ ; 
ELSE 
IF (SUM(K, X.L(I,K,T)) EQ 0, 
PUT “B-12°. 
ls 
es 


LOOP (K, 
TF (X-G(P ek Ty ro © 
PUT K.TL) ; 


Ve 


@15, LOOP(T, 
IF (ORD(T) GT 5, 
IF (ORD(T) LT 11, 
PUT TTL): 


LOGP(T,, .2Ur -/ 1 .2a; 
LOOP (T, 
IF (ORD(T) GT 5S, 
EF AORD (2) ET ii, 
IF (SHIPAVAL(I,T) EQ 0, 
PUT ‘UNAVAL’ ; 
ELSE 
IF (SUM(K, X.L(1I,K,T)) EQ O, 
POT 2 B24, 


LF (XL, KT soa, 
PUT K.TL)*; 


ie 
i 
lg 
PUT // 
@15, LOOP(T, 
IF (ORD(T) GT 10, 
PUT T.TL) ; 


LOOP PUL 7 1. Th: 
LOOP (T, 
IF (ORD(T) GT 10, 
IF (SHIPAVAL(I,T) EQ 0, 
PUT ’UNAVAL’ ; 
ELSE 
IF (SUM(K, X.L(1, 7) ) S500, 
PUT ‘’B-12’; 
3 
ae 


LOOP (K, 
IF (X.L(17 R79) bo 2 
PUT K.TL) ; 
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SONTEXT 

The following code creates the input file for interface with 
spread sheet. The name of the file is FAL1.PUT. 
SOFFTEXT 


FILE FAL1 ; 
PUT FAL1 ; 
FAL1.PC = 5 ; 
LOOP (T, 
LOOP (I, 
* LOOP (K, 
EE. (X2L (5) K, T) O01 
PUT o/s Ls Tie oe Pik. 2) 
} 
IF (SHIPAVAL(I,T) EQ 0, 
Pur / I-TL, TTL, “UNAVAL’ : 
ELSE 
IF (SUM (K, X.L(I,K,T)) EQ 0, 
PUT / I.TL, T.TL, ‘’B-12’; 
iF ° 


Ne 
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